The Social Behavior and the Evolution of Sexually Transmitted Diseases 
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We introduce a model for the evolution of sexually transmitted diseases, in which the 
social behavior is incorporated as a determinant factor for the further propagation of the 
infection. The system may be regarded as a society of agents where in principle anyone 
can sexually interact with any other one in the population. Different social behaviors are 
reflected in a distribution of sexual attitudes ranging from the more conservative to the more 
promiscuous. This is measured by what we call the promiscuity parameter. In terms of this 
parameter, we find a critical behavior for the evolution of the disease. There is a threshold 
below what the epidemic does not occur. We relate this critical value of the promiscuity 
to what epidemiologist call the basic reproductive number, connecting it with the other 
parameters of the model, namely the infectivity and the infective period in a quantitative 
way. We consider the possibility of subjects be grouped in couples. In this contribution only 
the homosexual case is analyzed. 

PACS numbers: 87.19.Xx, 87.23. Ge, 89.65.-s 



INTRODUCTION 



In the recent years, many mathematical models of so- 
cial phenomena have been formulated in order to describe 
a wide variety of phenomena 0. A particular interest 
was shown by the evolution of epidemic processes when 
the structure of the underlying society is taken into ac- 
count ||, m |], H . In many of these models the existence of 
a threshold value that determines the further evolution of 
a nucleus of infection was verified . On the contrary, 
it is possible to find other situations where the absence 
of threshold was shown |6). Nevertheless, the interest- 
ing feature in all of these models is that they incorporate 
a social aspect not taken into account in previous epi- 
demiological analysis. In the present work, we analyze a 
particular family of diseases, those sexually transmitted 
such as AIDS, Hepatitis B, Syphilis, etc. We will focus on 
a particular aspect of the society very related to the prop- 
agation of these diseases: the sexual behavior. We will 
consider a society composed by sexually active individu- 
als and with different behavioral patterns ranging from 
the most conservative or stable one, those who only have 
sexual intercourse with their stable and unique mate, to 
the most promiscuous who are continuously willing to 
change their sexual partner. The dynamic of the disease 
will be associated to the SIR case, where, at a given 
time, each individual in the population can be in one of 
the following three stages: susceptible (S), infective (J), 
and refractory or removed (R). A susceptible individual 
can become infective through contagion by an infective 
individual. Once an element has been infective, it en- 
ters a cycle that, after a fixed infection time, ends when 
the element reaches the refractory (or removed in case 
of death, as in the present model) state. At this stage 



the individual cannot be infected again nor infect the 
others, therefore the infection ultimately leads to defini- 
tive removal of elements from the susceptible population. 
This family of models has been used to describe the dy- 
namics of well-known infectious diseases, such as AIDS, 
rabies, and black death, by means of mean field approx- 
imation 01. 



THE MODEL 

Our model of the society conforms to a set of N sex- 
ually active subjects, which in principle are grouped in 
couples with the exception of a controlled proportion of 
singles (p sn g)- Initially a "promiscuity" value pi is ran- 
domly assigned to each agent from a semi-Gaussian dis- 
tribution of width p. pi determines the individual's ten- 
dency to dismiss his stable mate and go out -or just to 
go out in the case of the singles- to look for an occasional 
intercourse. In terms of the model the individual promis- 
cuity pi is the probability of trying to meet somebody else 
on each opportunity, determined by the time step of the 
simulation. Then, a susceptible individual can become 
infective, after a sexual intercourse with an infective one, 
with probability /3. The infective individuals remain in- 
fectious for a period r, after what they are removed by 
death. 

Every time step each individual must choose to have 
sex with his couple or with somebody else, the latter oc- 
curs with probability pi. Those who have decided to go 
out, must chose a partner at random and if the latter 
has decided to go out too, the occasional couple is made. 
Therefore there is no social structure a priori. The web 
of contacts is constructed dynamically during the simula- 
tion depending on the percentage of singles and p. Thus, 
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a random proportion of subjects, from those who have 
decided to go out, are reorganized in temporary couples 
every time step. Those who were abandoned by their 
couples and decided no to look for another partner will 
not have sex during this time step. The same is true for 
the singles who decided not to go out. Eventually, an in- 
dividual who goes out to look for somebody else may not 
have success remaining sexually inactive. In the present 
model, once an individual is removed, no new members 
are added to the society. This more complex situation 
will be analyzed in a further work. 

Summarizing, the model has four parameters: p sng , p, 
t, and (3 The first one controls the proportion of singles, 
the second one is the width of the semi-Gaussian distribu- 
tion of ^, while the latter's two are the same for all indi- 
viduals. As long as we assume that all agents are -except 
of the promiscuity and the marital status- identical (i.e. 
the infectivity between any two subjects is symmetrical), 
the model may be regarded as representing the spread 
of HIV/ AIDS in a male homosexual community with no 
distinction between passive and active partners. 

NUMERICAL RESULTS 

As was stated in the Introduction our results are for 
fixed number (N) of subjects, all sexually active -in prac- 
tice an individual could be sexually low active if single 
and with low pi. No new members are incorporated dur- 
ing the dynamics and no other cause of death apart from 
the sexually transmitted disease is considered. There- 
fore the number of susceptibles (S), infectives (I), and 
removeds (R) is conserved: N = S + I + R. Most of 
the simulation are performed for 100000 subjects popu- 
lation with only one individual initially infective. A hun- 
dred realizations (different runs starting with a different 
shuffling of the pi among individuals) were considered to 
make statistics. In order to give the simulations a human 
dimension one time step represents a day, and we make 
runs up to 5yr-100yr of simulation time, depending on 
the selected parameters (namely r). 

If we take for example p sng — (all subjects paired), 
p = 1/30, t = 1 yr, and (3—1 day^ 1 , the resulting time 
evolution of the infectives, I(t), and susceptibles, S(t), 
is similar to what can be obtained solving the classical 
SIR differential equations J^]. However, the similarities 
are only qualitative, because in our scheme different runs 
for the same parameters give different positions for the 
infectives maximum peak occurrence for example. Even 
more, a small percentage of the realizations does not pro- 
duce an epidemic; i.e the infection dies out taking just 
few individuals. This is an unexpected result from the 
point of view of mean field standard continuous models, 
however it does not seem unreasonable at all, consider- 
ing that people, behavior and randomness are involved 
in real life. Thus statistically there is a low but non zero 




promiscuity 



FIG. 1: Asymptotic percentage of death as a function of the 
promiscuity distribution width p. Parameters are r = lyr 
and /3 = 1. Notice the onset of the epidemic at the critical 
value p c = 0.9/30. Numerical results obtained with N = 1000, 
10000, and 100000 subjects, 100 realization and for p 3ng = 
case (all subjects paired). 

probability for the infection not to evolve with these pa- 
rameters, while the standard SIR model would give an 
epidemic as unique, deterministic outcome. Statistically, 
these cases should be considered as fluctuations. 

The results for several runs with the above parameters, 
but for p varying in the range 1/60 — 1/15, and for dif- 
ferent populations sizes are presented in Fig. [l| There it 
can be seen that the probability for an epidemic to occur, 
changes abruptly from zero to a finite value, at the critical 
promiscuity p c — 0.9/30. We can see too that there is a 
non-zero probability for non-epidemic at values of p > p c 
as denoted by the points at the horizontal axis, showing 
an asymptotic very low number of death. At the same 
time, the asymptotic number of death increases very fast 
with p > Pc, for realizations that produce epidemics. The 
line in Fig. [j] is the average of asymptotic deaths calcu- 
lated over all realizations which gave an epidemic (% of 
asymptotic death > 1%). In what follows we show this 
average instead of plotting the individuals runs. Proba- 
bly, the most interesting feature found in this work -not 
previously reported- is the epidemic threshold in terms 
of the promiscuity. The critical value of the population 
promiscuity, p c , depends obviously on the values of the 
other parameters (r, f3, and p sng ) as we will see. 

In the same Fig. |l| it can be noted that the obtained 
results are not an effect of the finite size of the sample. 
Several realizations on different sized population are plot- 
ted there, and the coincidence of all data is quite evident. 

Another aspect to take into account is the possibility of 
the epidemic evolution of being critically affected by the 
number of initially infected individuals, 1(0). Though it 
is clear that for a huge proportion of initially infectives, 
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the evolution will drastically change, allowing for a small 
increment in the number of them, namely from 1 to 5, in 
samples of 10000 subjects -and even in a small popula- 
tion of 1000 subjects- gives the same results. 

Up to here the results presented are for p sng = 0, i.e. 
the non singles case. How the results depend on this late 
parameter? Performing simulations for three different 
cases, i.e. p sng — 0, 0.5, 1, we can see that the qualitative 
behavior of the results is the same, with a displacement 
to higher values of p c as p sng increases. The results are 
plotted in Fig. ^, where the shift in the value of p c is 
observed, and also a slower evolution to high values of 
death as p increases. This is an expected result: once 
one of the members of the couple gets infective, sooner 
or later the other gets too; therefore, the more people are 
in couples, the more effective is the propagation of the 
epidemic. We have to remark here that in the present 
model the individual promiscuity is independent of the 
marital status. 
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FIG. 3: Squared critical value of the distribution promiscuity 
width (pc) as a function of the inverse product of f3 and r, 
for three different values of p ang ■ Promiscuity is a per month 
probability, where f3 is a per contact probability; r is in years. 



100 



in 

c3 

<U 

O 

o 



c/j 
c3 







1 


1 1 








/ 






1 ' ' 






1 / 








1 / 








/ / / 






/ ' ' 








f / 
















/ / 


— p . =0 








1 sing 






/ / 


-.p. =0.5 






' / 


r smg 






.-p. = 1 






1 1 


r sing 




1 L 


1 ; 


1 









1/30 2/30 
promiscuity 



3/30 



FIG. 2: Asymptotic percentage of death as a function of 
distribution promiscuity width for three different values of 
psng ((3 = 1 day -1 , r = 1 yr). Numerical results for 100000 
subjects and 100 realizations. 

It is known from classical SIR models (see ||) that 
the relevant parameter is p = a/r, where r is the infec- 
tiveness and a is the removal rate. We expect to verify 
the same for /3 and r, so the relevant parameter should 
be their product. This was numerically checked for dif- 
ferent values of (3 and r confirming the collapse of all 
data. Consequently the critical promiscuity depends on 
the product of (3 and r as can be observed in Fig ^, where 
p% as a function of (/3t) _1 is plotted along with linear fit- 
tings for the three values of p sng studied. From that we 
have obtained the expression p 2 c (3t = 0.36, 0.48, 0.65 (re- 
spectively for p sng = 0, 0.5, 1) connecting the p, (3, and r 
parameters at the critical point. 

As for practical purposes, we can test the threshold 
condition, obtained from our model, against HIV infec- 



tion data. Assuming an infectivity (3 — 0.1 |g] and a life 
time r = 5yr, we can roughly estimate p c — 0.045, this 
in turn means ((p) = (2/-k)^ 1 / 2 ^p for a semi-Gaussian dis- 
tribution) an average critical promiscuity of (p) = 0.036, 
i.e one "excursion" every 28 days in the population aver- 
age. Thus, this estimate is consistent with the outbreak 
of HIV infection among male homosexuals in the late 70's 
in the US. 

So far we have considered a semi-Gaussian distribution 
of promiscuity. It is known today that the web of sexual 
contacts has a power low distribution as shown in |l^] 
for a sample of the Swedish population, so the natural 
question is: Are all the results presented here valid for 
an actual society or are they mere academic speculation 
with no practical application? It has been put forward 
for example that the scale-free structure of the sexual net 
make it susceptible for an epidemic to occur without any 
threshold [|]. In order to verify this we have repeated 
some runs changing the semi-Gaussian distribution by a 
power-law distribution: 



n(p) 



if p < po 

2pyp~ 3 otherwise 



and we have found again a threshold for the epidemic 
outbreak For this distribution the parameter po has to 
be bigger than 0.1 for an epidemic to occur for (3 = 1 
and t = 1. The results can be observed in Fig. ||[ where 
we plot the asymptotic number of death for the power law 
promiscuity distribution together with the Gaussian dis- 
tribution. In order to compare both distribution we have 
normalized the promiscuity axis to the mean population 
promiscuity (p) . In this terms the onset of the epidemic 
are similar with a shift to smaller values in the power law 
case. Therefore the critical behavior is independent of the 
type of distribution. This can be explained because of the 
distribution cut-off; in the real world the net of contacts 
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have no nodes (subjects) with infinity connectivity. In 
our specific model, the maximum possible promiscuity is 
1 which means that the person goes out to meet someone 
everyday. Even in this extreme case he has a maximum 
for the number of possible encounters (less than 365 a 
year because he has to be accepted too). We think that, 
at the same time that our hypothesis it reasonably for 
a real population it gives an explanation for the appar- 
ent controversy between our results and the theoretical 
arguments against the existence of threshold 0] . 
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FIG. 4: Asymptotic percentage of death as a function of 
the mean promiscuity of the population (p) for a power-law 
distribution and a Gaussian distribution of promiscuity. Pa- 
rameters are r — lyr and (3 = \day~ 1 . Numerical results ob- 
tained with 100000 subjects, 100 realization and for p sng = 
case (all subjects paired). 



CONCLUSIONS 

We have presented a model for the spread of sexual 
transmitted diseases, taking into account the social be- 
havior. In spite of the simplistic way of modeling the 
social interactions, the model produce an outcome not 
previously reported in classical SIR models for infectious 
diseases of any kind, sexually transmitted or transmitted 
by other ways: a quantitative condition for the epidemic 
to occur. This threshold, called by epidemiologist, the 



basic reproductive number, was conjectured in the clas- 
sical literature in an empirical way Jl(], [ll]]. Moreover, it 
is not predicted, nor consistent with classical SIR mod- 
els. On the other hand the present model gives a direct 
connection between disease and social parameters and 
the epidemic threshold. Being quantitative this connec- 
tion could be of great utility to predict epidemics. The 
relevant parameters are four: two corresponding to epi- 
demiological aspects and related to the disease, the prob- 
ability of being infected in each sexual encounter (3 and 
the average life time of the infectious individual t, and 
two related to social aspects: the distribution of sexual 
contacts of the population pi, and the fraction of people 
that are not in stable couples. The knowledge of the last 
three would led us to determine the tolerated values for 
the infectivity /? to slow down and, why not, to stop the 
HIV infectious for example. 
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